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Abstract 
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1 Introduction 



It is well known that rare decays and/or oscillations of S-Mesons impose constraints on the 
parameter space of models Beyond the Standard Model (BSM): BSM contributions are not 
necessarily suppressed, once the dominant contributions both in the SM and BSM arise from 
loop diagrams (or are even absent in the SM). 

Recently, considerable progress has been made both on the experimental side (such as 
improved measurements of small branching ratios) and on the theoretical side, i.e. improved 
evaluations of SM predictions and BSM contributions. 

The purpose of the present paper is to study the resulting constraints on the parameter 
space of supersymmetric extensions of the standard model, both in the MSSM and the 
NMSSM, from BR{B ^^7), AM,, AM^, BR{Bs i^+i^-) and BR{B+ r+Ur). In 
the MSSM, similar analyses have recently been performed in [1-4] (see also refs. [5-10] for 
recent discussions within the Minimal Flavour Violating MSSM) . 

In [1, 2] the new experimental B physics results have been used to constrain the parameter 
space of the MSSM. In [3] it has been argued, that the new results on BR{B^ — > T'^Vr) are 
evidence for BSM contributions. A general fit has been performed in [4] in the context of 
the CMSSM (with universal Susy-breaking terms at the GUT scale) and the NUHM (with 
non-universal Higgs mass terms) , together with constraints on the dark matter relic density. 

One purpose of the present paper is to consider constraints from BR{B — > ^57) on 
the NMSSM. Our result is that the NMSSM specific effects on BR{B X^-f) are rather 
weak: in the NMSSM the charged Higgs mass squared receives (at tree level) a negative 
contribution relative to the MSSM which lowers its mass somewhat; once the result of 
BR{B X^-f) is plotted against Mh± , no difference between the MSSM and the NMSSM 
remains visible, however. Two loop corrections (relevant at large tan (3) are sensitive to the 
neutralino sector which includes the singhno in the NMSSM; we find, however, that even 
for relatively large singlino - MSSM-like-neutralino mixings the NMSSM specific numerical 
effect on BR{B — > Xg^) is numerically neghgible. (Combined constraints on the parameter 
space of the NMSSM from LEP, the dark matter rehc density and B physics - but without 
the recent developments in B physics - have been investigated previously in [11, 12].) 

Note that in the general MSSM, LEP constraints on the hghtest Higgs mass impose 
tan/3 > 3 (or tan/3 > 10 in the CMSSM). In the NMSSM (and the CNMSSM), LEP 
constraints on Higgs masses and couphngs allow for rather low values of tan /3 [13, 14]; here 
tan/3 can be as low as 1.5. Our results for BR{B — > Xg^) for low values of tan/3 (which 
have not been considered in [1-4]) are thus specific to the NMSSM, although the results in 
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the MSSM (without LEP constraints) would have been the same. 

In the NMSSM, important new contributions to B physics observables can originate 
from the presence of a relatively light CP odd Higgs boson [15-22], which can also be con- 
sistent with the dark matter relic density [11, 12], and which can contribute significantly via 
s-channel single and double penguin diagrams to B physics processes even for small tan/?. 
Constraints from B physics observables on this region of the parameter space of the NMSSM 
will be discussed in section 5. 

Our numerical results are obtained with the help of a Fortran code, that will be made 
public as a part of the NMSSMTools package [23]. It allows us to combine the constraints 
on the parameter space from B physics with constraints on the Higgs sector from LEP. (In 
the MSSM, subroutines that compute B physics observables are included in FeynHiggs [24], 
Suspect [25], MicrOmegas [26, 27] and Spheno [28]. Once all the calculations described below 
are included in NMSSMTools, it can also be used for the MSSM, since the MSSM is just a 
particular limiting case of the NMSSM.) 

In the remaining part of the introduction we briefiy review the experimental and the- 
oretical status of the various B physics observables, which are considered in the present 
paper. 

In the past constraints from b S'j have been particularly severe, since the experimental 
world average for BR{B ^s7) was somewhat below the (NLO) SM prediction [29,30], 
whereas at least the contribution involving a charged Higgs boson in the relevant diagram 
is positive. 

This situation has changed considerably during the last years: the present world average 
estimated by the Heavy Flavour Averaging Group [31] reads (for > Eq = 1.6 GeV) 

BRiB Xs-i) = (3.55 ± 0.241° ?^ ± 0.03) x 10"^ (1.1) 

exp 

The SM NNLO (0(a^)) corrections to the total BR{B Xs^t) branching fraction have 
recently be combined [32, 33], which give 

BR{B ^^7)1^^ = (3.15 ± 0.23) x 10"^ (1.2) 

In [34] the treatment of the cut E^ > 1.6 GeV on the photon energy has been improved, 
leading to a still lower SM prediction: 

BR{B ^s7)|^^^ = (2.98 ±0.26) x 10"^ (1.3) 

This result can be interpreted as (still weak) evidence for BSM contributions to 6 — > S7; in 
any case constraints on the parameter space of Susy models have become considerably less 
stringent. 
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Next we turn to ^Mg^d- AM^ has recently been measured by the CDF collaboration [35] 
with the result 

AMf = 17.77 ± .12 ps-^ . (1.4) 

A standard model prediction 

AMf ^ = 20.5 ± 3.1 ps-^ (1.5) 

can be obtained using a determination of |VjsVtf,| = (41.3 ± .7) x 10~^ from tree level pro- 
cesses (where effects from BSM physics affect the higher order corrections only) [6], and a 
determination of /s,-/^ = 0.281 ± .021 GcV by the HPQCD collaboration [36]. (In [36], 
the central value AMf^ = 20.3 ps"^ has been obtained, since iV^j^V^bl = 41.0 x 10~^ has 
been used. We note that here and below the CKM matrix elements are defined in terms of a 
low energy effective Lagrangian, whose parameters are determined from low energy processes 
[37]. In [37], these CKM matrix elements are denoted by but we omit the subscript 
" e//" in the following.) Hence, a negative contribution to AMg from BSM processes would 
be welcome. 

AMd is quite well known [31], 

AMf ^ = 0.507 ± .004 ps'^ . (1.6) 

Again, a standard model prediction (see also [8]) 

AMf^ = 0.59 ±0.19 ps-^ (1.7) 

can be obtained using a determination of |V^^l^t6| = (8. 6 ±1.4) x 10^'^ from tree level processes 
[6], as above, and /b./^/Zb^v^ = 1.216 ± .041 from [38]. 

The various Susy diagrams which contribute to AMg (g = s, d) are box diagrams involving 
charged Higgs bosons, stops and charginos (see, e.g., [39]), and double penguin diagrams 
involving neutral CP even or CP odd Higgs bosons whose contributions increase like tan^ j3 
for large tan/? (see [37] for a detailed analysis). As a function of the mass M^j of the Higgs 
boson, these contributions to the Wilson coefficients behave like 1/Mfj, and depend on the 
mixing angles of the CP even and CP odd Higgs mass matrices. In the MSSM, the dominant 
contributions ~ ^/Mfi (where h denotes the lightest Higgs scalar) cancel at large tan/? [37], 
and one is left with contributions ~ 1/^1 (where A denotes the CP odd scalar in the MSSM, 
whose mass is close to the heavy CP even scalar for large Ma) which cannot be too large, 
given the lower bound on in the MSSM. 

In the NMSSM, three neutral CP even and two CP odd Higgs bosons (we neglect the 
Goldstone boson here) contribute to the double penguin diagrams. Notably the lightest 
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CP odd Higgs boson Ai can be quite light in the NMSSM and escape the present LEP 
constraints [15-21], but with couphngs strong enough to generate large effects for low Ma-^ 
[17]. Interestingly, the resulting contributions to AMg are negative which can improve the 
agreement with its measurement. 

Neutral Higgs bosons with effective ffavour violating couplings contribute also to 
BR{Bs iJ,'^ii~), where the new CDF result is at 95% confidence level [40] 

BR(Bs a-^a-) < 5.8 x 10"^ . (1.8) 

exp 

(At present, constraints from BR{Bd — > jJ-'^jJ-') are less restrictive.) The SM prediction is 
still smaller by an order of magnitude [41,42], 

BR{Bs ^ /^V") = (3.8 ± 0.1) X 10-^ , (1.9) 

which leaves some room for BSM contributions. Again, a light CP odd Higgs boson Ai can 
lead to an important effect in the NMSSM [17]; in the case of BR{Bs jl^jjT), however, 
its contribution must not be too large. 

Finally we turn to BR{B'^ — > r+i/,-), which has been observed by the Belle [43] and 
BABAR [44] experiments. The actual world average performed by the Heavy Flavor Aver- 
aging Group [31] is 

BR(B+ T+i/^) = (1.32 ± .49) x 10"^ . (1.10) 

exp 

Unfortunately, the corresponding SM prediction is handicapped by a large uncertainty 
concerning the CKM matrix element |\45| [31,45]: Its determination from inclusive semi- 
leptonic b decays gives values near \Vub\ ~ 4.4 x 10^'\ whereas its determination from 
exclusive semileptonic decays gives values near \ Vuh\ ~ 3.7 x 10^'^ (leading to a discrepancy of 
the order of 2a). Accordingly, together with the uncertainties from the hadronic parameter 
Jb, quite different SM predictions for BR{B^ t^^t) can be obtained, ranging from 
BR{B+^T+Ur) = {0.85 ±.13) xlQ-'^ [1] to BR{B+ T+Ur) = (1.59 ± .40) x 10"^ 



SM 



SM 



[3]. Hence we will allow for quite large theoretical error bars on this process with the result 
that it hardly constrains the Susy parameter space; this situation can change in the future, 
however. 

In section 2 we describe the sources of the contributions to BR{B ^s7) that we 
take into account. In section 3 we give the sources of our calculations of AM^, AM^, 
BR{Bs fJ'^fJ' ) and BR{B^ r'^u^). In section 4 we present results for BR{B X^'-f) 
both for large tan f3 (relevant for the MSSM and the NMSSM) and for low tan /5 (relevant 
for the NMSSM only). 
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Finally, in section 5, we investigate combined constraints from BR{B Xs^), AMg^d, 
BR[Bs fJ'~^fJ' ) and BR{B^ t^Vt) in parameter regions relevant simultaneously for the 
MSSM and the NMSSM, and also on the NMSSM specific region involving a light CP odd 
Higgs scalar. In all cases we include constraints on the parameter space from LEP on Higgs 
masses and couplings as in the updated version of NMHDECAY [46,47]. In section 6 we 
conclude with a summary and an outlook. 



2 Computation of BR{B X,7) 



The starting point of our computation is the expression for the branching ratio as in [29, 30], 

vt,*' subtracted 



BR (b Xsj 
BR (b X^eu) 



exp 



Vcb 



6ae 



\K, + r{iio)Kt + e,^f + B{Eo) + N{Eo) , (2.1) 



valid for a matching scale /iq — mt{mt)^^ . 
In (2.1), we use [33] 

BR {b Xceu 

and 



exp 



0.1061 , 



ch 



= 0.967 



from tree level processes [6] . 
C in (2.1) is given by 
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for which we use the numerical value 



C = 0.580 . 



(2.2) 
(2.3) 



(2.4) 



(2.5) 



Eq is the lower cutoff on the photon energy, for which we chose Eq = 1.6 GeV. 
includes the SM top quark and the BSM contributions, whereas Kc denotes the SM charm 
quark contribution. r(/Xo) is the ratio m^'^(/xo)/m^^, for which we use [29, 30] 

-0.003 



r{iio) = 0.578 



0.1185 / V4.69 



1.25 



/ /^o 
U65 



-0.08 / \ 0-006 

469 



(2.6) 



with ml^ = 4.68 GeV as in [33] and /ib = mb{mb). (The dependence on the scale is in 
fact negligibly small.) 
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In (2.1) denotes the electroweak radiative corrections, B[Eq) the (gluon) brems- 
strahlung corrections, and N{Eq) are nonperturbative corrections. 

Strictly speaking, the expression (2.1) is vahd to NLO, where the charm quark contribu- 
tion {Kc) can be separated from the top quark/BSM contribution {Kt). depends on the 
ratio nic/mb, and hence on the scheme and the scale at which these masses are taken. On 
the one hand this ambiguity is a NNLO effect, which is responsible for the largest part of 
the theoretical error in the NLO result [30] 

BR{B ^ X,7) "^^"^ = (3.6ll°i^ ^ ^ ± .02cKM ± 0.24^„. ± 0.14,ea/e) X 10-\ (2.7) 

mc/mb 

We found that the NNLO result (1.2) is reproduced (for rrit^poie — 171.4 GeV, as assumed 
in [32, 33]), if one uses the relatively large value 

777 

— = 0.307 (2.8) 
nib 

(close to the pole quark masses) in the expression for Kc- We beheve that as long as the 
BSM contributions - which are added hnearly to the SM contributions in the factor Kt - are 
not evaluated to NNLO, the error arising from this procedure is not larger than the error 
intrinsic to the BSM contributions (which is estimated quite conservatively below). It is 
guaranteed, in any case, that the result for the BR{B — > ^57) in the decoupling hmit of the 
BSM contributions assumes the NNLO SM value in (1.2). 

Subsequently we describe the origin of the formulas used for our evaluation of the quan- 
tities Kc, N{Eq), B{Eo), eew and Kt in (2.1). First, Kc is computed as in Eq. (3.7) in [29], 
with /lb — TTib, the value (2.8) for nic/nib and 

/xo = mtinit)^ (2.9) 

for the matching scale /iq. The ratio of CKM matrix elements e^, that appears in Eq. (3.7) 
in [29], is taken from [33]: 

e, = CK6/(V^:i/t6) = -0.011 + i 0.0180 (2.10) 

The nonperturbative corrections N{Eq) are computed as in Eq. (3.10) in [29] in terms of 
the lowest order coefficients K^^ and Kf^^ (including the BSM contributions to the latter), 
with A2 = 0.12 GeV^. {N{Eq) is actually independent from Eq in this approximation). 

The bremsstrahlung corrections B{Eq) are taken from the appendix E in [29] with, we 
repeat, an energy cutoff Eq = 1.6 GeV. For the parameter z = (mc/mbY we use a value 
consistent with eq. (2.8) above. (In any case the dependence of B{Eq) on z is weak [29].) 
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The corrections ~ (with q = s) as in Eq. (28) in [30] are taken into account, with given 
in Eq. (2.10) above. The contributions to B{Eq) from the coefficients Cf^ with i = 3. . .6 
are neglected as in [29], on the other hand the BSM contributions to the coefficients C'^^ 
and Cf^ are taken into account. 

For the electroweak corrections e^w in (2.1) we use the formula (3.9) in [29] (see also Eq. 
(27) in [30]), which gives a SM contribution ef^ = 0.0071 according to [48]. To this SM 
value for we add the BSM contributions as in [29, 30] in terms of the BSM contributions 
to the coefficients Cy^g discussed below. 

Finally we turn to the calculation of Kt including the BSM contributions. First, the SM 
contributions to Kt (including the NLO in ag) are taken from Eq. (3.8) in [29], with the 
above Eq. (2.9) for the matching scale /xq- The BSM contributions are added as in Eq. 
(5.1) in [29]. The BSM contributions appear in the LO Wilson coefficients C'^^^^'^ {jio), 
^{o)BSM ^^^^ and the NLO Wilson coefficients C^^^^^^ {/j^q), C^^^^^^^/J^o) and C^^^^^^ {/j^q) of 
the corresponding operators Pj. 

Our calculation of these Wilson coefficients within the MSSM and the NMSSM starts 
with the calculation of the corrections e^, e[ and e[ to the couplings of the charged Higgs 
bosons to quarks defined in [49] (see also [50]), which are important at large tan/?. We use 
the expressions for these parameters given in [27], which include sbottom and electroweak 
contributions, and in which a sign error in [49] is corrected. In the case of and e[ we sum 
over the 5 neutralino states of the NMSSM with its corresponding masses and couplings. 
(In the MSSM limit A, k ^ of the NMSSM, the fifth neutralino decouples and does not 
contribute.) Then we proceed with the computation of the following BSM contributions to 
the Wilson coefficients: 

a) The chargino-squark loop contributions to Cj^^ and Cg°^ (as, e.g., in appendix B in 
[27]), computed again at Mgusy and evolved to our matching scale /io. Corresponding NLO 
corrections are known in the particular case where one stop is lighter than the other squarks 
and the gluino [51], and the complete QCD corrections have been computed in [52], but here 
we content ourselves with the summation of the leading logarithms of the ratio Msusy/ l^o 
via the RG evolution of the Wilson coefficients. 

b) The charged Higgs-top-quark loop contributions to C'f^ and C'^'' (again as in appendix 
B in [27]), and the corresponding NLO contributions to C^l\ C^"* and Cg^"* [53]. The LO 
contributions to Cf^ and Cf^ are evolved from the scale corresponding to the charged Higgs 
mass to our matching scale /^o, and we took care not to include large logarithms - that appear 
potentially also in the NLO contributions - twice. (Higher order large tan j3 corrections to 
the NLO contributions are neglected.) 
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c) As in [54] we take the neutral Higgs contributions to the Wilson coefficients Cy and 
Cg^"* into account following Eq. (6.61) in [37]. However, contrary to AMg and BR{Bs — ^ 

below, these neutral Higgs effects remain small and usually inside our theoretical 
error bars. 

d) Finally the large tan j3 corrections induce also a shift in the SM contributions to the 
coefficients Cf'\jj,Q) and Cf\iJ,o) [27,53]. 

Herewith we have described completely the origins of the considered contributions to 
BR{B ^ X,7). 

3 AM^, BR{Bs and BR{B+ r+z/^) 

In this section, we discuss the sources for our evaluation of the B physics observables AM^ 
(g = s, d), BR{Bs — > A*"*"/^") and BR{B^ — > t'^i't)- The formula for AMq is taken from [37], 
eqs. (6.6-7): 

F^t\ (3.1) 
with 

F^, = So{x,) + ^Cll:^ + " (Cf ^ + Cf^«) + P^^C^^ + . . . (3.2) 

where we have omitted neglicibly small contributions, and where we take [37] r = 0.985, 
p5LL ^ _Q pLR ^ Q 9Q ^ Q 55;^ ^Yie meson masses Mb^ = 5.2794 GeV 

and Mb,, = 5.3675 GeV, and the hadronic parameters /b^ \J Bb, = 0.281 GeV from [36] and 
\l Bb^ = 0.231 GeV from Jb, \J Bb,/ \l Bb^ = 1-216 [38]. As stated in the introduction, 
we use the CKM factors deduced from tree level processes, which are less sensitive to BSM 
physics: |V;*,\46| = 41.3 x lO'^ and \V*^Vtb\ = 8.6 x lO^^ [6]. 5*0 in Fu stands for the SM 

contribution [xt = ( 1 ), whereas the coefficients contain BSM contributions to the 
corresponding effective 4-quark operators. 

Let us discuss the various contributions to F/J which we take into account (we repeat that 
we assume minimal flavor violation such that the only source of flavor violation is the CKM 
matrix): The SM contribution originates from quark/W^^ box diagrams. In multi-Higgs 
extensions of the SM such as the MSSM or the NMSSM, charged Higgs bosons can replace 
one or both bosons in these box diagrams. A second type of box diagrams arises in Susy 
from squark/chargino loops. All these box contributions are calculated as in eqs. (93-95) in 
[39] and added directly to 5*0: 

So^So + Xt{AH±+Al^) (3.3) 



Mb^VbILBb, 
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We have checked that at low tan /3, where the box contributions are most significant, the 
resuhs in [10] are reproduced. 

Double Penguin diagrams involving a neutral Higgs propagator connecting two flavor 
changing effective vertices can be significantly enhanced for large tan j3 or light scalars. We 
closely follow the analysis carried out in [37]: 

• First, we compute flavor dependent e parameters (effective vertices) arising from loops 
involving sparticles in the effective Lagrangian describing the Higgs quark couplings. 
We use Eq. (5.1) and appendix A. 2 in [37]. However, we extend the neutralino sector 
according to the NMSSM; the corresponding generalization of the MSSM formulae is 
straightforward. 

• Next, we deflne flavor- changing neutral Higgs-quark couplings ^f^/^j^ as in Eqs. (3.55- 
56) in [37] {S denote the various neutral Higgs bosons). The corresponding Higgs 
mixing angles and can be generalized in a straightforward way to the NMSSM 
using the decomposition of the neutral weak eigenstates and into the neutral 
physical states (in the convention of [37]) as H^* = J2so — J2s° ^d^^- 

• Finally, we use Eq. (6.12) of [37] (neglecting the Goldstone boson contribution) for 
the three relevant coefficients Cf^^, Cf^^ and €2^. However, as we will face very 
hght (pseudo)scalar masses (possibly below 10 GeV in some parts of the NMSSM 
parameter space), we can no longer be content with the approximation ^ for the 
scalar propagator (see [17], Eq. (32)). Thus, we replace these factors by Breit-Wigner 



(The width Ts is computed as in NMSSMTools [46,47].) In the MSSM, relations between 
the Higgs masses at large tan f3 allow for further simplifications of the final formula for AM,^ 
(see [37], Eq. (6.23)). However, in the NMSSM a correct description of the singlet like 
contributions does not allow for such simplifications. 

Next we consider BR{Bs — > ix^ix~). We calculate the Branching Ratio according to Eq. 
(5.15-16) of [55] (we neglect the cQ: 



functions: 




sgn{ml - Ml) 



(3.4) 






CA (3.5) 



2 
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where ca contains the SM contribution arising from box and penguin diagrams, which is one 
order of magnitude below the sensitivity of present experimental data. The neutral Higgs 
contributions to cs and cp are the only ones which could lead to a significant deviation from 
the SM prediction. The corresponding diagrams involve the effective ffavour violating neutral 
Higgs vertex and a neutral Higgs propagator. We infer from an appropriate generalization 
of the equations given in [37] the appropriate formulae for the coefficients cs and cp in the 
NMSSM. Again, it proves necessary to replace the approximation ^ by a Breit-Wigner 
function. 

Charged Higgs corrections to BR{B^ T'^z^r) were studied in [56] and lead to a de- 
structive interference with the SM (H^"'') contribution: 

BR{B^ r+..) = (l-^Xfs K^frnVH , (3.6) 



where th parametrizes the deviation from the SM prediction. The expression for th has 
been improved in [57] in order to take large tan/3 corrections into account: 



1 - 



\mH±J l + iotanP 



(3.7) 



Having described the origin of all relevant calculations, we turn to the numerical results, 
concentrating first on BR{B — > ^^7). 



4 Results for BR{B X^j) in the MSSM and the 
NMSSM 

The BSM contributions to B ^ Xgj depend essentially on the charged Higgs mass, tan /? 
and, for large tan/3, on Af. 

First, we focus on the impact of the charged Higgs mass on BR{B — > X57), which is 
always positive. The branching ratio is a decreasing function of mH±, since the contribu- 
tions from charged Higgs diagrams decay fike l/m^±. Before the recent improvements on 
the experimental side and the SM contributions discussed in the introduction, quite severe 
bounds on m/f± could be deduced notably for small to modest values of tan/3, where the 
additional Susy contributions (which can have both signs, depending on the relative sign of 
At to ijl) cannot be too large in absolute value. 

The updated situation is described in Figs. 1-4. In Fig. 1 we show our results for 
B — > Xg^ for tan/3 = 5, universal squark masses of 1 TeV, gaugino masses Mi = 150 GeV, 
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Figure 1: BR{B — > Xg'-f) as a function of the charged Higgs mass, for tan/3 = 5, At = 
±2500 GeV. The green hues represent the experimental la bounds. 

Ms = 300 GeV, M3 = 900 GeV, for two extreme values of At = +2.5 TeV and -2.5 TeV as 
a function of m//±. We scan over the parameter /i between +100 GeV and +1 TeV, which 
explains the broadening of the two dotted distributions. (The inner regions correspond 
to larger values of /i, the outer regions to the lowest value of fi that is allowed by the non- 
observation of charginos.) For the top quark mass we take 171.4 GeV. The la experimentally 
allowed region is also indicated and it becomes clear that, at least after taking theoretical 
errors into account (see below), relatively low values of mH± down to ~ 200 GeV cannot be 
excluded. This result holds both for the MSSM and the NMSSM (where the //-parameter 
has to be replaced by an effective parameter /Xg// = A {S), we use the conventions of [46]); 
no dependence on the additional parameters of the NMSSM remains visible. 

Before we turn to larger values of tan/5, we study a relatively low value tan/? = 2.2 
which would make it very difficult for the MSSM to satisfy the constraints from LEP on 
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Figure 2: BR{B — > Xg'-f) as a function of the charged Higgs mass, for tan/5 = 2.2, At = 
±2500 GeV. The green hues represent the experimental la bounds. 

the hghtest neutral Higgs boson mass, but which is perfectly consistent in the NMSSM [13], 
even in the CNMSSM with universal soft terms at the GUT scale [14]. Fig. 2 is the same as 
Fig. 1, but for tan/3 = 2.2 and NMSSM parameters A = 0.5, k = 0.4 and = -200 GeV, 
which lead to neutral Higgs masses consistent with LEP bounds provided > 300 GeV 
(due to correlations between the various Higgs mass matrices in the NMSSM). There is no 
particular impact of the NMSSM parameters on B ^sl, however. One finds that this 
NMSSM specific region in parameter space is hardly constrained by this observable. 

Next we investigate B — >■ Xg'y for larger values of tan (3. We find an approximate linear 
dependence on tan j3 with a slope determined essentially by At, at least for given /i, which we 
fix now at 300 GeV. In Fig. 3 we show our results for various values of At, = 300 GeV 
(and the same other parameters as above), and in Fig. 4 for m//± = 1 TeV (which is obtained 
essentially by a vertical shift of Fig. 3). Now one finds that, the larger tan/5, the stronger 
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Figure 3: BR{B Xg'-f) as a function of tan (3, for Mj^h 
At 



300 GeV and various values of 



are constraints on l^^l from BR{B Xg'-f). On the other hand, positive contributions from 
relatively light charged Higgses can easily be cancelled by appropriate contributions from 
squark/chargino loops. Again, these results hold both for the MSSM and the NMSSM. We 
note that all points/lines shown in our Figures correspond to parameters which satisfy LEP 
constraints on Susy Higgs bosons, but this is not always trivial: small values of 1^4^! and 
tan/? can lead to a too light neutral Higgs boson both in the MSSM and in the NMSSM; 
this is the reason why we confined ourselves to \At\ > 600 GeV in Figs. 3 and 4, and why 
the lines (notably for \At \ = 600 GeV) do not continue to arbitrarily small values of tan/5. 

Although further dependencies on, e.g., the soft Susy breaking squark and gaugino masses 
would certainly merit further studies (which can be performed using NMSSMTools [23], once 
updated), we believe that our Figs. 1-4 represent a fairly comprehensive review of the actual 
status of the predictions for BR{B -> Xs'j) in the MSSM and the NMSSM. 
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Figure 4: BR{B —>■ Xg'-f) as a function of tan (3, for Mh+ = 1 TeV and various values of At 



5 Constraints from BR{B -^s7), BR{Bs fi^ii ), 

AMg, and BR{B+ t^p^) in the MSSM and the 
NMSSM 



The aim of this chapter is to study the combination and the relative relevance of the 
constraints on the parameter space of the MSSM and the NMSSM from BR{B ^sl), 
BR{Bs —>■ fi'^fi^), AMg^d and BR{B^ t^Ut-). To this end we need to estimate the theo- 
retical error implicit in our calculations. We intend to remain conservative and to denote a 
point as excluded only if one of the observables falls outside the 95% confidence limit (or 2a). 

In the case of BR{B — > Xg'j), the theoretical error will depend on the parameters of 
the Susy model under consideration; a general value for the theoretical error would be 
misleading. Hence we estimate the theoretical error separately for the charged Higgs, Susy 
and SM contributions to BR{B X^'j) as follows: Since the charged Higgs contribution 
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is evaluated to NLO, we assume that its relative theoretical error is only 10%. For the 
Susy contribution, which is evaluated to LO only (up to leading logarithms), we assume a 
(conservative) relative theoretical error of 30%. Finally we estimate the theoretical error 
bars of the SM contribution as follows: Given that the improved treatment of the cut on 
the photon energy in [34] leads to a lower SM prediction than in [32], we allow the SM 
contribution to BR{B -> Xs^f) to vary in the range 2.72 x lO"'' to 3.38 x 10""^. The SM and 
BSM errors are added linearly, which gives our estimate of the final theoretical error. 

For AMg and BR{Bs fi~^fi~), we estimate the theoretical error due to BSM contribu- 
tions to be of the order of 30%, since no QCD corrections are taken into account. We add 
these uncertainties linearly to the 2a SM error bars, which gives our complete theoretical 
error estimate. (The Icr SM error bars on AMg, arising mostly from the uncertainties of 
CKM matrix elements and lattice computations of hadronic parameters, are given in eqs. 
(1.5) and (1.7) above.) 

Concerning BR{B~^ t~^i>t-), the uncertainties originating from the CKM matrix ele- 
ment \Vub\ are considerable. We allow \Vub\ to vary in the range 3.3 x 10~^ < \Vub\ ^ 4.7 x 
10~^, with 4.0 X 10"^ as central value. For Jb we use Jb = 0.216 ± 0.022 GeV as obtained 
by the HPQCD collaboration [58]. It just so happens that the corresponding central values 
lead to a SM prediction BR{B~^ — > t'^'^t) = 1-32 x 10~^ in agreement with the experi- 
mental central value given in (1.10). Allowing for 2a error bars on and the experimental 
average (1.10), and using the above range for \Vub\ one finds that th in Eq. (3.7), neglecting 
additional theoretical errors, is allowed to vary over the quite large range 

0.13 < r^^ < 4.0 . (5.1) 

Consequently the constraints on the parameters tan (5 and mj^± from this process are typi- 
cally less stringent than the ones from other processes. 

Now we turn to the dependence of the observables on the most relevant parameters. 

AMq and BR{Bs fJ'^fJ' ) are quite sensitive to (double) Penguin contributions involv- 

ing neutral Higgs bosons. These contributions are controlled by the parameter I — j , 
where z/ = 3 for BR{Bs — > A^^Ai") and u = 2 for AMg (and denotes the component 
of the neutral Higgs boson S); this explains why BR{Bs fJ'^fJ' ) is usually more sensitive 
to neutral Higgs effects, at least at large tan/? where they can become huge, leading to a 
violation of experimental bounds both in the MSSM and the NMSSM. 

Thus, in general, large values of tan/5 are rather strongly constrained by these observ- 
ables. However, it is still possible to reduce the neutral Higgs contributions by assuming 
heavy scalars and pseudoscalars (through a large doublet mass Ma ~ Mh±)- Another way to 
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circumvent these constraints consists in assuming parameters as the trihnear soft-couphng 
At or fieff such that the e parameters (which control the flavour violating neutral Higgs 
couplings) remain small enough - here cancellations are often possible. 

Only for low tan l3 can the positive contributions from Susy box diagrams to AMg be more 
important than the double Penguin contributions. Given that the SM prediction for AM^ 
[36] is already ~ la above the CDF result [35], such additional positive BSM contributions 
could in principle exclude points in the parameter space at low tan f3. (For larger tan (3 
the double Penguin diagram, which gives a negative contribution AMg, usually dominates 
the box diagrams.) However, once we use 2a error bars for the CKM matrix element and 
hadronic uncertainties, such exclusions at low tan (3 occur scarcely in practice. 

In the following we present several examples of constraints on tan (3 and Mh± (for fixed 
other parameters) that originate from the above processes. 

First we consider the MSSM and the NMSSM for relatively small values of A and k, ( < 
0.1), for which the contributions to the above processes are practically the same in both 
models. For the soft Susy breaking squark and gaugino masses we take the same values as 
in Figs. 1-4, and 300 GeV for /i (or fi^jf in the NMSSM). 

In Fig. 5 we assume At = 2.5 TeV. Dark dotted regions are excluded by LEP: Here 
and in Figs. 6 and 7 below the non-observation of a light neutral Higgs scalar h at LEP 
implies lower limits on the MSSM parameter M4 (depending on tan/? and At) which, in 
turn, lead to lower limits on Mh± (~ Ma for large Ma)- The domain allowed by LEP 
is further constrained by B physics processes. We note that the BR{B Xg'j) is by far 
the most stringent constraint in Fig. 5. It is indeed particularly severe because both the 
charged Higgs and the Susy contributions are positive and thus cannot balance each other. 
(Constraints from AMd are never more restrictive than constraints from AMg, hence AMq 
in Figs. 5 and 6 stands for AMg.) 

In Fig. 6 we switch to At = —2.5 TeV, where the situation is quite different (the notation 
is the same as in Fig. 5): BR{B Xg'^) allows for additional domains, which originate from 
cancellations between the charged Higgs and Susy contributions (strongly enhanced by the 
large value of \At\). Therefore, light charged Higgs bosons (with masses down to ~ 100 GeV) 
are not excluded by this process; on the contrary, for tan j3 > 20, they must be light enough 
to avoid a large decrease of the branching ratio due to Susy diagrams. However, these 
regions are also constrained by BR{Bg — > n^n^) and, less stringently, by BR{B^ — > t^i't) 
and AMq. 

In Fig. 7 we consider smaller values of At = 600 GeV: Now, small values of tan f3 and 
Mh± (or Ma) are ruled out by LEP constraints on Mh. (The precise bound is very sensitive 
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Figure 5: Constraints in the tan (3-Mh+ plane for At = 2500 GeV. 

to radiative corrections to and hence to rutop- We recall that we use = 171.4 GeV.) 
LEP constraints do not rule out a narrow strip around Mh± ~ 120 GeV (already visible 
in Fig. 6), where the coupling of h to the Z-Boson is suppressed (since the MSSM-like 
parameter sin(/5 — a) happens to be small) and where Mh ~ 100 GeV. However, even this 
region is now excluded by the charged Higgs contribution to BR{B —>■ Xg'j). (For positive 
or small absolute values of At the Susy contribution to BR{B — > Xs'~f) cannot cancel the 
charged Higgs contribution.) BR{Bs —>■ /U^yU~) does no longer lead to constraints since 
neutral Higgs effects, which are (roughly) proportional to At, remain small for a low value 
of this parameter. On the contrary, BR{B^ t^Ut), which depends only weakly on At, 
can become the dominant B physics constraint for tan (3 > 30. 

Next we discuss a region specific to the NMSSM: In the NMSSM, singlet-like pseu- 
doscalars Ai even below 10 GeV are able to survive LEP constraints. However, their loop 
induced flavour violating couplings to quarks and leptons can be large enough to cause sig- 
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Figure 6: Constraints in the tan (3-Mh+ plane for At = —2500 GeV. 

nificant contributions to B physics observables, most of all for CP odd scalar masses near 
the resonance (m^^ ~ ^b,) and/or large tan/5. (Now, large values of tan/5 do not only lead 
to larger couplings of the light CP odd scalars, but also to an increase of their width which, 
in turn, enhances their contribution via the s-channel Penguin diagram even for masses a 
few GeV away from the resonance.) 

In the following we present several examples of constraints on tan j3 and M^^ that orig- 
inate from B physics processes. For the soft Susy breaking squark and gaugino masses we 
take the same values as above, and 300 GeV for /ie//- The NMSSM specific parameters are 
chosen as A = 0.45, k = 0.4 and = —30 GeV. However, Ax (or the MSSM-like parameter 
M\ = XS {Ax + kS) / {cos P sin (3)) must be chosen within a relatively narrow tan /^-dependent 
window such that LEP constraints on all CP even and CP odd Higgs scalars remain satisfied. 
In Figs. 8-10 is chosen within this ~l-2 GeV wide window. {Ma varies from 300 to 
400 GeV for tan/5 ~ 1.5 to 10; LEP constraints would allow to extend this window up 
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Figure 7: Constraints in the tan j3-MH+ plane for At = 600 GeV. 



to tan (3 = 50 with more finetuning on M^', however, B physics constraints exclude this 
domain.) 

In Fig. 8 we consider the plane Ma^ vs. tan/3, and assume At = —2.5 TeV. Now, the 
constraints from BR{Bs —>■ /U^/U") are the most relevant, and lead to strong upper limits 
on tan/5 at fixed M^^. Among the remaining observables, constraints from AMg are also 
significant but generally redundant with the respect to BR{Bs /i"*"/!"). 

Whereas the situation for At = 2.5 TeV in Fig. 9 is similar to the one with At = —2.5 TeV 
(the main difference comes from BR{B —* Xg'y), which excludes now a region with very light 
Ai and tan/3 > 7 already covered by BR{Bs — > ji^jj,^)), the case At = 500 GeV considered 
in Fig. 10 is quite different: In contrast to Figs. 8 and 9 the lightest scalar Higgs mass rrih 
is below 90 GeV, but LEP constraints can still be satisfied due to the decay h —* AiAi. 

Note that, on the one hand, Ai in Fig. 10 has a ~ 90% singlet component, but also 
a ~ 40% doublet component. For tan/3 near 5 its coupling to down type quarks is even 
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Figure 8: Constraints in the tan /J-Ma^ plane for At = —2500 GeV. 

~ 2 times larger than the one of a SM scalar Higgs boson. As function of Ma^ (and rrih), 
LEP constraints on h ^ AiAi ^ 4 b, h ^ AiAi ^ A t or h —* AiAi 4 jets have then 
to be applied, which explains the jumps in the upper bound on tan f3. However, within the 
region allowed by LEP, B physics constraints are particularly weak: only a narrow stripe 
with near M^^ is excluded by BR{Bs fJ'^fJ' ) and AMg. Once again, this is due to 
the fact that neutral Higgs effects are essentially proportional to \At \ and small for small \At\. 

6 Summary and Outlook 

In this article, we have updated constraints from B physics observables on the parameters 
of the MSSM and the NMSSM (assuming minimal flavour violation), combining them with 
LEP constraints on the parameter space. Available SM and BSM radiative corrections are 
included in the calculations, which will be made public in the form of a Fortran code. 
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Figure 9: Constraints in the tan/9-M^j plane for At = 2500 GeV. 

As expected, constraints from BR{B —>■ Xg'j) have become weaker due to the recent 
increase of the world average, and the decrease of the SM prediction (which is now below 
the experimental average). Our numerical results (summarized in Figs. 1-4) show that 
constraints still arise if, simultaneously, Mh± is small {Mh± < 300 GeV) and tan (3 not too 
large ( < 10), or if tan/? > 10 and \At\ is large. We have verified explicitely (for the first 
time), that NMSSM specific contributions to BR{B Xg'y) are numerically negligible. 

Among the other processes, BR[Bs —>■ fi~^fj^') is typically the most sensitive and can 
exclude regions in parameter space for tan/? > 15 that would be allowed by BR{B Xg'y), 
see Fig. 6. However, also BR{B^ t^Vt) can lead to the most relevant constraints for very 
large tan /5, cf. Fig. 7. 

In the NMSSM specific case of a light CP odd Higgs scalar, constraints from BR{Bs 
fj'~^fj'~) (inside the LEP allowed region) are quite strong for large \At\ (cf. Figs. 8 and 9), 
but exclude only a small region around Ma^ ~ 5 GeV for small \At\ (cf. Fig. 10). 
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Figure 10: Constraints in the tan plane for At = 500 GeV. 

In the future, our calculations will allow to combine constraints from B physics observ- 
ables with additional assumptions such as universal soft terms at the GUT scale (the CMSSM 
and the CNMSSM) and/or constraints from the dark matter relic density via NMSSM- 
Tools [23]. 
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